Quantum entanglement produced in the formation of a black hole 
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A field in tlie vacuum state, which is in principle separable, can evolve to an entangled state in a 
dynamical gravitational collapse. We will study, quantify, and discuss the origin of this entanglement, 
showing that it could even reach the maximal entanglement limit for low frequencies or very small 
black holes, with consequences in micro-black hole formation and the final stages of evaporating 
black holes. This entanglement provides quantum information resources between the modes in the 
asymptotic future (thermal Hawking radiation) and those which fall to the event horizon. We will 
also show that fermions are more sensitive than bosons to this quantum entanglement generation. 
This fact could be helpful in finding experimental evidence of the genuine quantum Hawking effect 
in analog models. 
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I. INTRODUCTION 

Quantum entanglement has been recognized to play a 
key role in black hole thermodynamics and the fate of 
information in the presence of horizons; some previous 
studies were performed in stationary cases, namely the 
eternal acceleration scenario and the stationary eternal 
Schwarzschild black hole [i|-Qj ^ot addressing issues re- 
lated with dynamics and time evolution of gravitating 
quantum fields. On the other hand some studies involv- 
ing entanglement dynamics in expanding universe sce- 
narios have shown that the interaction with the gravita- 
tional field can produce entanglement between quantum 
field modes [i,[l^. 

In this paper we analyze the issue of entanglement pro- 
duction in a dynamical gravitational collapse. With this 
aim, we consider both a bosonic (scalar) and a fermionic 
(Grassmann scalar) field which initially are in the vac- 
uum state and compute their asymptotic time evolution 
under the gravitational interaction in a stellar collapse. 
The vacuum state evolves to an entangled state of modes 
in the future null infinity (which give rise to Hawking 
radiation [ll|) and modes that do not reach it because 
they fall into the forming event horizon. 

We will argue that the initial vacuum state in the 
asymptotic past does not have any physical quantum en- 
tanglement, and that it evolves to a state that is phys- 
ically entangled as a consequence of the creation of the 
event horizon. This entanglement depends on the mass 
of the black hole and the frequency of the field modes. In 
particular, for very small frequencies or very small black 
holes, a maximally entangled state could be produced. 

The entanglement generated in a gravitational collapse 
thus appears as a quantum resource for non-demolition 
methods aiming to extract information about the field 
modes which fall into the horizon from the outgoing 
Hawking radiation. These methods would be most rele- 
vant for cases such as the formation of micro-black holes 



and the final stages of an evaporating black hole when 
the mass is getting smaller and, therefore, quantum cor- 
relations generated between the Hawking radiation and 
the infalling modes grow to become even maximal, as we 
will show. 

Earlier works proved that fermions and bosons have 
qualitatively different behaviors in phenomena such as 
the Unruh entanglement degradation Q, 0, 0] and the 
entanglement generation in the background of expand- 
ing universes [9|, llOj . Here, we will show that for fermions 
the generation of entanglement due to gravitational col- 
lapse is more robust than for bosons. This robustness is 
more evident from the peak of the thermal spectrum of 
Hawking radiation towards the ultraviolet. 

Previous works in the literature (see for example |12l - 
|l5| among many others) showed that Hawking radiation 
is correlated with the field state falling into the collaps- 
ing star. However neither the analysis of the associ- 
ated entanglement entropy as a function of the black 
hole parameters nor the comparison between fermionic 
and bosonic behavior have been carried out so far. The 
study of these issues, the nature of the entanglement 
produced in a gravitational collapse and, more impor- 
tant, its dependence on the nature of the quantum field 
(bosonic/fermionic) is decisive in order to gain a deeper 
understanding about quantum entanglement in general 
relativistic scenarios as it was proven for other setups 
such as acceleration horizons, eternal black holes and ex- 
panding universes [3, 0, 0, 0, [l3 • 

Since entanglement is a pure quantum effect, under- 
standing its behavior in these scenarios can well be rele- 
vant to discern the genuine quantum Hawking radiation 
from classical induced emission in black hole analogs [ia| 
(see, for example, Ref. [l3), where both classical and 
quantum perturbations obey the same evolution laws. 
It will also follow from our study that fermionic modes 
could be more suitable for this task since they are more 
reliable in encoding entanglement information. 





FIG. 1: Carter-Penrose diagrams for gravitational collapse: 
Stellar collapse (left) and Vaidya spacetinie (right). 

Finally, we will argue that the entanglement between 
the infalling and the Hawking radiation modes neither ex- 
isted as a quantum information resource nor could have 
been acknowledged by any observer before the collapse 
occurs, namely in the asymptotic past. This is impor- 
tant in order to understand the dynamics of the creation 
of correlations in the gravitational collapse scenario since 
these correlations are exclusively due to quantum entan- 
glement, as discussed in the literature [l3 - [l5| . 



II. GRAVITATIONAL COLLAPSE 

In order to analyze the entanglement production in- 
duced by gravitational collapse we will consider the 
Vaidya dynamical solution to Einstein equations (see e.g. 
Rcf. [13,]) that, despite its simplicity, contains all the in- 
gredients relevant to our study. Refinements of the model 
to make it more realistic only introduce subleading cor- 
rections. The Vaidya spacetime (Fig. [1]), is conveniently 
described in terms of ingoing Eddington-Finkelstein co- 
ordinates by the metric 

ds^ = -{l-^MM\dv^ + 2dvdr + r^dn^, (1) 

where r is the radial coordinate, v is the ingoing null co- 
ordinate, and M{v) ~ mO[v — vq). For vq < v this is 
nothing but the ingoing Eddington-Finkelstein represen- 
tation for the Schwarzschild metric whereas for v < vq 
it is just Minkowski spacetime. This metric represents a 
radial ingoing collapsing Shockwave of radiation. As it 
can be seen in Fig. [H Uh ~ vq — Am represents the last 
null ray that can escape to the future null infinity J^^ 
and hence that will eventually form the event horizon. 

Let us now introduce two different bases of solutions 
to the Klein-Gordon equation in this collapsing space- 
time. On the one hand, we shall define the 'in' Fock 
basis in terms of ingoing positive frequency modes, asso- 
ciated with the natural time parameter v at the null past 
infinity J^^, which is a Cauchy surface: 
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On the other hand, we can also consider an alternative 
Cauchy surface in the future to define another basis. In 
this case, the asymptotic future ^+ is not enough and we 
also need the future event horizon ^+. The 'out' modes 
defined as being outgoing positive-frequency in terms of 
the natural time parameter T^out at ^+ are 
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where r^out = v — 2r*^t and r*^^ is the radial tortoise 
coordinate in the Schwarzschild region. It can be shown 
(see e.g. Ref. [l^) that, for late times ryout — J- oo at J^+, 
these modes u°"' are concentrated near v-a at ,f^ and 
have the following behavior: 
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These modes have only support in the region v < Vn- 
This is evident as only the light rays that depart from 
V < Vii will reach the asymptotic region ^+ since the 
rest will fall into the forming horizon defined by u = Uh. 
This is the only relevant regime, as far as entanglement 
production is concerned. 

For the 'hor' modes defined at J^^, there is no such 
natural time parameter. A simple way to choose these 
modes is defining them as the modes that in the asymp- 
totic past J^~ behave in the same way as m°"* but de- 
fined for V > Vy,, that is to say, as modes that leave the 
asymptotic past but do not reach the asymptotic future 
since they will fall into the horizon. This criterion is the 
simplest that clearly shows the generation of quantum 
entanglement between the field in the horizon and the 
asymptotic region. In any case, since we will trace over 
all modes at the horizon, the choice of such modes does 
not affect the result. Therefore, we define the incoming 
modes crossing the horizon by reversing the signs oivji — v 
and oj in (j4]) so that, near J^~ these modes are 
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We are now ready to write the annihilation operators 
of bosonic field modes in the asymptotic past in terms 
of the corresponding creation and annihilation operators 
defined in terms of modes in the future: 
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where tanhr^^ = g^'*'^™". The precise values of ^p and 
a^^i are not relevant for this analysis. 

Hence the vacuum 10);^^, annihilated by ([6]) for all fre- 
quencies uj' , acquires the following form in terixis of the 
'out-hor' basis: 



|0)i„ = iVexp j^tanhr^ 




where A^ — (ricj <^o^h''tj) is a normalization constant. 
We can rewrite this state in terms of modes \n^) with 
frequency a; and occupation number n as 
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III. ANALYZING ENTANGLEMENT 

This is a two-mode squeezed state. Therefore, it is 
a pure entangled state of the modes in the asymptotic 
future and the modes falling across the event horizon. 
Given the tensor product structure no entanglement is 
created between different frequency modes. Hence, we 
will concentrate the analysis in one single arbitrary fre- 
quency w. 

We can compute the entropy of entanglement for this 
state which is the ultimate entanglement measure for a 
bipartite pure state, defined as the Von Neumann entropy 
of the reduced state obtained upon tracing over one of the 
subsystems of the bipartite state. To compute it we need 
the partial state pout = trhor(|0)i„(0|), which turns out to 
be Pout = Ylu Pout,u, where 
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This is, indeed, a thermal radiation state whose temper- 
ature is nothing but the Hawking temperature (87rm)~^, 
as it can be easily seen. However, this is only the par- 
tial state of the field, not the complete quantum state, 
which is globally entangled. If we compute the entropy 
of entanglement Se,cj = tr(pout,Ld log2 Pout,^) for each fre- 
quency, after some calculations, we obtain 



Se,lj = (coshr^^) log2 (coshr^) 
- (sinhr,^)^ log2 (sinhr^) , 



(10) 



which is displayed in Fig. [2] As (|S]) is a pure state, all the 
correlations between modes at the horizon and modes in 
the asymptotic region are due to quantum entanglement. 
Analogously we can compute the entanglement for 
fermionic fields. If we consider a spinless Dirac field (ei- 
ther one dimensional or a Grassmann scalar), the anal- 
ysis is entirely analogous considering now both particle 
and antiparticle modes. We assume again that the ini- 
tial state of the field is the vacuum that, after some long 
nontrivial calculations, can be expressed it in the Fock 
basis at the asymptotic future and the 'hor' modes: 

|0)in = n[(cO«''-)'|00)hor|00)o,t 
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FIG. 2: Entanglement between bosonic (continuous blue) and 
fermionic (red dashed) field modes in -j^^ and in y'^. The 
lesser the mass of the star or the mode frequency, the higher 
the entanglement reached. 



where tanfj^ = g-4irmw^ Here, we are using the double 
Fock basis, the first figure inside each ket representing 
particles and the second antiparticles. 

We can compute the entropy of entanglement of this 
pure state. The partial density matrix in the asymptotic 
future Pout = trhor(|0)i„(0|) = H^ Pont.u,-, is given by 

;; ^4 
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S \2 



(sin2f^) 



|1^0)„J1^0| + |01J„J01, 



-|-(sinf;^) |lc^li^)out(l"l"l 



(12) 



which is again a thermal state with Hawking temperature 
(Sttto)"^, and 

5E,a; = -2[(cosf„)^log2(cosf„)^+(sinf„)^log2(sinf^)^], 

(13) 
which is also displayed in Fig. [51 

Figure [2] shows that the entanglement decreases as the 
mass of the black hole or the frequency of the mode 
increase. When comparing bosons with fermions one 
must have in mind that the entropy of entanglement is 
bounded by (the logarithm of) the dimension of the par- 
tial Hilbert space ('out' Fock space in our case). There- 
fore, due to Pauli exclusion principle, the maximum en- 
tropy of entanglement for fermions is Se,u: ~ 2, which cor- 
responds to a maximally entangled state. On the other 
hand, for bosons, the entanglement is distributed among 
the superposition of all the occupation numbers and the 
entropy can grow unboundedly, reaching the maximally 
entangled state in the limit of infinite entropy. In this 
sense, the entanglement generated in the fermionic case 
is more useful and robust due to the limited dimension 
of the Fock space for each fermionic mode. 

This result can be traced back to the inherent differ- 
ences between fermions and bosons. Specifically, it is 
Pauli exclusion principle which makes fermionic entangle- 
ment more reliable. Similar results about realiability of 
entanglement for fermions were also found in the expand- 
ing universe scenarios [10| . This responds to the high in- 
fluence of statistics in entanglement behaviour in general 



rclativistic settings as it was investigated in [y, Q • On the 
other hand, Vaidya space-time has all the fundamental 
features of a stellar collapse and shows how the entangle- 
ment is created by the appearance of an event horizon. 
Hence, in other collapsing scenarios or including the sub- 
leading grey-body factor corrections, these fundamental 
statistical differences will not disappear. The qualita- 
tively different behavior of entanglement for bosons and 
fermions is not an artifact of choosing a particular col- 
lapse scenario but is due to fundamental statistical prin- 
ciples. 

In the above analysis we have considered plane wave 
modes, which are completely delocalized. However, an 
entirely analogous analysis can be easily carried out using 
very well localized Gaussian states, with the same results 
about quantum entanglement behavior. 



IV. CONCLUSIONS 

We have shown that the formation of an event hori- 
zon generates entanglement. If we start from the vac- 
uum state in the asymptotic past, after the gravitational 
collapse process is complete we end up with a state in 
the asymptotic future which shares pure quantum cor- 
relations with the field modes which fall into the hori- 
zon. One could think that this entanglement was already 
present before the collapse, arguing that (as proved in 
[iSl p the vacuum state of a quantum field can be un- 
derstood as an entangled state of space-like separated 
regions. In other words, if we artificially divided the 
Cauchy surface in which the vacuum state is determined 
into two parts, we would have a quantum correlated state 
between the two partitions. In principle we could have 
done a bipartition of the vacuum state in J^~ such that 
it would reflect entanglement between the partial state 
of the vacuum for w < Uh and the corresponding partial 
state for v > Vh. However, it is not until the collapse 
occurs that we have the information about what Va is. 
So, achieving beforehand the right bipartition (trying to 
argue that the entanglement was already in the vacuum 
state) would require a complete knowledge of the whole 
future and, consequently, there is no reason 'a priori' to 
do such bipartition. The entanglement, eventually gen- 
erated by the collapse, will remain unnoticed to early 
observers, who are deprived of any means to acknowl- 
edge and use it for quantum information tasks. It is 
well known that if we introduce artificial bipartitions of 
a quantum system, its description can show entangle- 
ment as a consequence of the partition. However, not 
being associated with a physical bipartition this entan- 
glement does not codify any physical information and, 
hence, cannot be used to perform any quantum infor- 
mation processes. (One example of this kind of non- 
useful entanglement is statistical entanglement between 
two undistinguishable fermions [l9|). 

Gravitational collapse selects a specific partition of the 
initial vacuum state by means of the creation of an event 



horizon. In the asymptotic past there was no reason 
to consider a specific bipartition of the vacuum state, 
whereas in the future there is a clear physically meaning- 
ful bipartition: What in .y~ was expressed as a separable 
state, now becomes expressed in terms of modes that cor- 
respond to the future null infinity and the ones which fall 
across the event horizon. This means that gravitational 
collapse defines a particular physical way to break the 
arbitrariness of bipartitioning the vacuum into different 
subsystems. This gravitational production of entangle- 
ment would be a physical realization of the potentiality 
of the vacuum state to be an entangled state and is there- 
fore a genuine entanglement creation process. 

We have computed the explicit functional form of this 
entanglement and its dependence on the mass of the black 
hole (which determines the surface gravity). For more 
complicated scenarios (with charge or angular momen- 
tum), it will depend on these parameters as well. 

For small black holes, the outgoing Hawking radiation 
tends to be maximally entangled with the state of the 
field falling into the horizon for both bosons and fermions. 
This means that if a hypothetical high energy process 
generates a micro-black hole, a projective measurement 
carried out on the emitted radiation (as, for instance, the 
detection of Hawking radiation) will 'collapse' the quan- 
tum state of the field that is falling into the event horizon 
and give us certainty about the outcome of possible mea- 
surements carried out in the vicinity of the horizon. Fur- 
thermore, at least theoretically speaking, the available 
quantum information resources would be maximum and, 
therefore, one could perform quantum information tasks 
such as quantum teleportation with maximum fidelity 
from the infalling modes to the modes in the asymp- 
totic future y^ if the observer of the infalling modes 
managed to dispatch an outgoing classical signal before 
crossing the horizon. On the other hand, low frequency 
modes become more entangled than the higher ones. So, 
the infrared part of the Hawking spectrum would provide 
more information about the state at the horizon than the 
ultraviolet. 

Arguably, similar conclusions can be drawn for the fi- 
nal stages of an evaporating black hole: As the mass 
of the black hole diminishes, the temperature of the 
Hawking radiation spectrum increases, and therefore, the 
quantum state of the field tends to a maximally entangled 
one in the limit of m — > 0. 

We have seen that the entanglement generated in 
fcrmionic fields is more robust than for bosons. Although 
the entropy of entanglement in the zero mass limit is 
greater in the bosonic case due to the higher dimension 
of the partial Hilbert space, we have argued that the in- 
formation is more reliably encoded in the limited Fock 
space of fcrmionic fields. Furthermore, as we consider 
higher frequency modes, fermionic entanglement proves 
to be much more easily created by the collapse. What is 
more, the turning point in which the entropy of entan- 
glement for fermions becomes numerically larger than for 
bosons is actually near the peak of the thermal emission 



(Fig. ^. This means, that, in general, a measurement 
carried out on Hawking radiation of fermionic particles 
will give us more information about the near-horizon field 
state. This might also be useful in analog gravity real- 
izations as we have already discussed, specifically in sys- 
tems where the field excitations are fermionic (see e.g. 
Ref. [201), which would be, as shown, at an advantage 
over the bosonic cases. To account for this quantum en- 
tanglement in analog experiments one should carry out 
measurements of the quantum correlations between the 
emitted thermal spectrum and the infalling modes and 
detect Bell inequalities violations. This is easier as it 
gets closer to the maximally entangled case. 
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